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ABSTRACT 

A microscopic formulation of Haldane's exclusions statistics is given in terms 
of a priori occupation probabilities of states. It is shown that negative probabilities 
are always necessary to reproduce fractional statistics. Based on this formulation, 
a path-integral realization for systems with exclusion statistics is derived. This 
has the advantage of being generalizable to interacting systems, and can be used 
as the starting point for further generalizations of statistics. As a byproduct, the 
vanishing of the heat capacity at zero temperature for exclusion statistics systems 
is proved. 
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Statistics is an inherently quantum mechanical property of identical particles 
which, as the name suggests, modifies the statistical mechanical properties of large 
collections of such particles. It enters through symmetry properties of the wave- 
functions of many-body states, and this was the starting point, until some time 
ago, for various generalizations away from the standard bosonic and fermionic 
cases, such as parastatistics and anyons. Haldane, however, taking the term statis- 
tics more literally, defined a generalized exclusion statistics through the reduction 
of the Hilbert space of additional particles in a system due to the ones already 
present in the system [1]. He proposed then the definition 



where N is the number of particles in the system, d is the dimensionality of the 
single-particle Hilbert space, obtained by holding the quantum numbers of N — 1 
particles fixed, and Ad and AiV are their variation keeping the size and boundary 
conditions of the system fixed, g = corresponds to bosons (no exclusion) while 
g — 1 corresponds to fermions, excluding a single state for the remaining particles, 
the one they occupy. 

On the basis of (1), Haldane proposed the combinatorial formula for the number 
of many-body states of N particles occupying a group of K states 



Based on the above, Wu derived the thermodynamical properties of particles with 
exclusion statistics [2] (see also [3,4]), and this system has received a lot of re- 
cent attention [5,6,7]. Appropriate generalizations for several species of identical 
particles also exist. 

It is obvious from (1) and (2) that exclusions statistics makes sense only in a 
statistical sence, since Ad and M can become fractional for AiV = 1 or TV = 1. It 
is, nevertheless, useful to attempt a microscopic realization and interpretation of 
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fractional exclusion statistics, and see what it implies for the one-particle states, 
and this is what will be done in this paper. Such a description has the obvious 
advantage of being generalizable to interacting particles, for which the notion of d 
becomes hard to define. 

The starting point will be the grand partition function for exclusion-statistics 
particles ("g-ons," as called in [6]) in K states 



where we put x = exp(/i — e)/kT with \x the chemical potential and assumed that 
all K states are at the same energy e. In the statistical limit of large K, Z(K) 
should be extensive. This intruduces, then, the notion of a microscopic description 
of the system in which the above Z is the K-th power of a single-state partition 
function. Each single level can be occupied by any number of particles, but with 
an a priori probability P n for each occupancy n independent of the temperature. 
We thus demand 



for all x. The above probabilities must, in general, depend on K in order to satisfy 
(4). This reflects the fact that g-ons are not well-defined for microscopic systems 
(see also the remarks in [6]). If, however, P n (K) assume some (finite) asymptotic 
values as K goes to infinity (as they should for an extensive Z), then the above 
microscopic partition function becomes an accurate description in the statistical 
limit. To calculate P n in this limit, we first notice that the combinatorial formula 
(2) counts (at least for integer g) the ways of placing N identical particles in K sites 
arranged into a one-dimensional open lattice, under the restriction that any two 
particles be at least g sites apart. Clearly, for large K, the "boundary condition" 
that the lattice is open cannot influence the statistical mechanics of the system. 
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We choose, then, to examine instead particles placed on a periodic lattice under 
the same restriction. This modifies the combinatorics into 

_ K[K-(g-l)N-l]\ 
M ~ N\{K-gN)\ {b) 

Clearly (5) reproduces the standard bosonic and fermionic results for g = 0,1. 
Repeating the analysis of [2], it can be verified that M' leads indeed to the same 
statistical mechanics as M. It is now pointed out that the P n defined in terms of 
M' are independent of K. (Proof: Putting P n (K) = P n + 0(1/ K), where P n are 
the asymptotic values at K — oo, and equating terms of ra-th order in x we obtain 

KP n (K) + K(K - ^Pn-^P^K) + ■■■ = M'(K, n) (6) 

Assuming that all P m (K) for m < n are independent of K, all the terms in 
(6) other than KP n (K) are polynomials in K without constant term, and so is 
M'(K,n). Therefore P n (K) cannot contain any terms of 0(1/ K). Since P\ = 1, 
we inductively showed that all P n are i^-independent.) Therefore, P n can be 
calculated from M'(l,n), and we obtain 



Pn = 

It can be a postoriori checked that the expressions for P n (K) obtained from M 
converge to (7) for K = oo. The above P n for n — 0, . . . , 5 and g — \ agree with 
the values calculated in [6] using a different approach. 

The most obvious feature of the above expressions is that, unless g — 0, 1, they 
always become negative for some values of n. Therefore, their interpretation as 
probabilities is problematic. This is an inherent problem of fractional g-on statistics 
which cannot be rectified by, e.g., truncating M(K,N) to zero for iV > K/g. 
The description of the statistical system in terms of effective negative microscopic 
probabilities is, nevertheless, accurate and useful. Note, also, that the above P n 
never truncate to zero for n above some maximal value (unless g — 1), unlike 
parafermions. 
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From the above, the single-level partition function Z(x) = Z can be shown to 
satisfy 

z g _ z g-i = x ( 8 ) 

This is a transcedental equation which in principle determines Z and whose power- 
series solution reproduces the P n as coefficients. The average occupation number 
n is expressed as 

n = ^xd x Z = xd x W (9) 

where W = \nZ is the free energy (over —kT). It can be shown that (9) together 
with (8) imply for n 

(1 - gn) 9 [l - (g - ^n] 1 ' 9 = nx' 1 (10) 

in accordance with the result of [2,3,4]. 

One immediate consequence of the above relations is the vanishing of the zero- 
temperature heat capacity of a g-on system Co, defined as 

oo 

Co = j{e- V)dp[n(f3) + n(-(3) - -] (11) 



where f3 = 1/kT, and 1/g is the saturation density for n at zero temperature and 
e < /i. Using (9) we can express Cq as 

C = ^(a; = A)-ilnA-^(a; = 0) (12) 

where A is a cutoff to be taken to infinity. From (8) we can deduce that W(x = 0) = 
and W(x = A) = ± In A + O(A-Vfl). Therefore C = for all g, as conjectured 
in [6]. This expresses the fact that the ground state of the many-body g-on system 
is nondegenerate. This is expected as a generic feature of particle systems, but is 
explicitly verified here for g-ons. 
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It is easy to derive a duality relation for Z: 

Z-\g,x-n + Z-\- g ,x) = l (13) 

From the above relation and (8), (9), the duality relation for the density is recovered 
[4,6] 

gn(g, x) + -n(-, x^ 1 ' 9 ) = 1 (14) 

We regard the formula (13) as more fundamental since it seems to be more generic. 
For instance, parafermions of order p = 1/g are defined such that at most p particles 
can be put per state with probabilities 1. Thus 

1 _ x p+i 

Z par =l + x+---X P = — (15) 

1 — X 

from which we can write the generalized parafermionic partition function Z par (g, x) 
by simply putting p = 1/g above. It can be seen that Z par also satisfies (13), but 
not (14). 

The free energy W can be expressed as a power series in x 

oo 

W = J2—x n (16) 

n=l 

in terms of the "connected" weights 

wi = PI , w 2 = 2P 2 - Pi , w 3 = 3P 3 - 3PiP 2 + Pi (17) 
etc. We find for w n : 

m=l 

These are remarkably similar to P n (except for the range of m). Notice that the 
w n are not probabilities, but rather virial coefficients. In fact, w n — 1 for bosons 
and w n = (— ) n_1 for fermions. Also, w 2 = 1 — 2g [5]. 
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From the above expressions for w n we can find a path integral representation 
for the partition function of g-ons in an arbitrary external potential. We start 
from the usual euclidean path integral with periodic time (5 for iV particles with 
action the sum of iV one-particle actions, and sum over all particle numbers N with 
appropriate chemical potential weights. Since the particles are identical, we must 
also sum over paths where particles have exchanged final positions, with weights 
equal to the inverse symmetry factors of the parmutation to avoid overcounting 
(compare with Feynman diagrams). Thus the path integral for each N decomposes 
into sectors labeled by the elements of the permutation group Perm(N). By 
the usual argument, the free energy will be given by the sum of all connected 
path integrals. It is obvious that these are the ones where the final positions of 
the particles are a cyclic permutation of the original ones (since these are the 
only elements of Perm(N) that cannot be written as a product of commuting 
elements). These have a symmetry factor oil/N corresponding to cyclic relabelings 
of particle coordinates (compare with the factors of 1/n included in (16)). They 
really correspond to one particle wrapping N times around euclidean time. Thus, 
if we weight these configurations with the extra factors w^, as we have the right 
to do since they belong to topologically distinct sectors, we will reproduce the free 
energy of a distribution of g-ons on the energy levels of the one-body problem, that 
is 

oo „ N 

WGM = E e " N N / WN II Dx n (t n )e- S ^^ (19) 

N=l n=l 

where Se is the one-particle euclidean action and the paths obey the boundary 
conditions x n {(5) = x n+ i(0), xn(/3) = xq(0). (x can be in arbitrary dimensions.) 
The partition function will be the path integral over all disconnected components, 
with appropriate symmetry factors and a factor of w n for each disconnected n- 
particle component. 

It is clear that the above path integral is not unitary, since the weights w n are 
not phases, nor does it respect cluster decomposition, since the w n do not provide 
true representations of the permutation group (unlike the g — 0, 1 cases). This is 
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again a manifestation of the non-microscopic nature of exclusion statistics. It does 
make sense, nevertheless, at the statistical limit. 

The above path- integral realization can be extended to other statistics. E.g., 
for parafermions of order p (with p integer) the corresponding weigts w n are 

w n = — p for n = modfjo + 1) , 1 otherwise. (20) 

This representation is more economical than the one calling for p distinct flavors 
of fermions and projecting all states transforming in an irreducible representation 
of Perm(p) into a unique quantum state. The origin of the apparent non-unitarity 
and breakdown of cluster decomposition in the above integral for parafermions is 
clear: it is due to the above projection, which must be inserted in the (unitary) 
many-flavor path integral. 

The possibility to define statistics throught the choice of the coefficients w n 
suggests other possible generalizations. Perhaps the simplest one is to choose 

w n = (-a)"" 1 (21) 

that is, one factor of —a for each unavoidable particle crossing. This leads to the 
statistical distribution for the average occupation number n 

n = -i \a ( 22 ) 

which is the simplest imaginable generalization of the Fermi and Bose distribution. 
The combinatorial formula for putting N particles in K states for the above a- 
statistics is 

M = a N ( ^ )! = K(K - a)(K - 2a)---(K - (N - l)a) 

N\(£-N)\ N\ 1 1 

This can be thought as a different realization of the exclusion statistics idea: the 
first particle put in the system has K states to choose, the next has K — a due to 
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the presence of the previous one an so on, and dividing by N\ avoids overcounting. 
Fermions and bosons correspond to a — 1 and a = — 1 respectively, while a = 
corresponds to Boltzmann statistics (as is also clear from the path integral, in 
which no configurations where particles have exchanged positions are allowed, but 
factors of 1/N\ are still included). The corresponding single-level probabilities are 

n— 1 

^ tt 1 — mot , ,n 

Pn = 11 TT — 24 
A \ l + m 

m=l 

For a — 1/p with p integer (a fraction of a fermion), the above probabilities are 
all positive for n up to p and vanish beyond that. For a < all probabilities 
are positive and nonzero. Thus, the above system has a bosonic (a < 0) and 
a fermionic (a > 0) sector, with Boltzmann statistics as the separator. It is a 
plausible alternative definition of exclusion statistics, due to (23), and has many 
appealing features, not shared by the standard (Haldane) exclusion statistics, such 
as positive probabilities, a maximum single-level occupancy in accordance with the 
fraction of a fermion that a represents, and analytic expressions for all thermody- 
namic quantities. It would be interesting to find a physical system in which these 
statistics are realized. 

We conclude by pointing out that, once we have the path integral (19) we can 
easily extend the notion of exclusion statistics to interacting particles: we simply 
replace the action Yln^E[ x n\ by the full interacting iV-particle action, thus cir- 
cumventing all difficulties with combinatorial formulae. In the interacting case one 
has to work with the full partition function, rather than the free energy (19), since 
topologically disconnected diagrams are still dynamically connected through the 
interactions and do not factorize. Applications of the above on physical systems, 
as well as possible generalizations to the many-flavor mutual-statistics case are left 
for future work. 

Acknowledgements: I would like to thank D. Karabali and P. Nair for interesting 
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